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Abstract 



G : 

3 , 

cm : 

t— I , A phenomenological analysis of lifetimes of bottom and charmed hadrons 

, within the framework of the heavy quark expansion is performed. The 

baryon matrix element is evaluated using the bag model and the nonrela- 
\Q . tivistic quark model. We find that bottom-baryon lifetimes follow the pat- 

tern r(f2b) ~ T C^b) > T (A-b) — r (^b)- However, neither the lifetime ratio 
r(Af,) /r{Bd) nor the absolute decay rates of the baryon and B mesons 
, can be explained. One way of solving both difficulties is to allow the pres- 

, ence of linear 1 jvtiQ corrections by scaling the inclusive nonleptonic width with 

Q^h, the fifth power of the hadron mass tuhq rather than the heavy quark mass 

Qh. m Q- The hierarchy of bottom baryon lifetimes is dramatically modified to 

r(Aft) > r(E b ) > r(S^) > r(Q&): The longestdived Qb among bottom baryons 
^> [ in the OPE prescription now becomes shortestdived. The replacement of tuq by 

TriH Q in nonleptonic widths is natural and justified in the PQCD-based factor- 
ization approach formulated in terms of hadron-level kinematics. For inclusive 
charmed baryon decays, we argue that since the heavy quark expansion does 
not converge, local duality cannot be tested in this case. We show that while 
the ansatz of substituting the heavy quark mass by the hadron mass provides 
a much better description of the charmed-baryon lifetime ratios, it appears 
unnatural and unpredictive for describing the absolute inclusive decay rates of 
charmed baryons, contrary to the bottom case. 



PACS numbers: 12.39.Hg, 12.39.Jh, 13.20.He, 13.30-a 



1 



I. INTRODUCTION 



The lifetime differences among the charmed mesons D + , D° and charmed baryons have 
been studied extensively both experimentally and theoretically since late seventies. It was 
realized very early that the naive parton model gives the same lifetimes for all heavy par- 
ticles containing a heavy quark Q and that the underlying mechanism for the decay width 
differences and the lifetime hierarchy of heavy hadrons comes mainly from the nonspectator 
effects like W-exchange and Pauli interference due to the identical quarks produced in heavy 
quark decay and in the wavefunction (for a review, see [1,2]). The nonspectator effects were 
expressed in eighties in terms of local four-quark operators by relating the total widths to the 
imaginary part of certain forward scattering amplitudes [3-5]. (The nonspectator effects for 
charmed baryons were first studied in [6].) With the advent of heavy quark effective theory 
(HQET), it was recognized in early nineties that nonperturbative corrections to the parton 
picture can be systematically expanded in powers of 1/mg [7,8]. Subsequently, it was demon- 
strated that this 1/rriQ expansion is applicable not only to global quantities such as lifetimes, 
but also to local quantities, e.g. the lepton spectrum in the semileptonic decays of heavy 
hadrons [9]. Therefore, the above-mentioned phenomenological work in eighties acquired a 
firm theoretical footing in nineties, namely the heavy quark expansion (HQE), which is a 
generalization of the operator product expansion (OPE) in 1/ijiq. Within this QCD-based 
framework, some phenomenological assumptions can be turned into some coherent and quan- 
titative statements and nonperturbative effects can be systematically studied. As an example, 
consider the baryon matrix element of the two-quark operator (A b \bb\A b ). The conventional 
quark- model evaluation of this matrix element is model-dependent: 

(AfrNAfe) = | 1 NQM; 
2m Ab \ J d 3 r [ujj(r) — v\{r)\ bag model, 

where u(r) and v(r) are the large and small components, respectively, of the quark wavefunc- 
tion. However, the matrix element (1.1), which is equal to unity in the nonrelativistic quark 
model (NQM), becomes smaller in the bag model due to the contribution from the lower 
component of the quark wavefunction. In the HQE approach, it is given by [see Eq. (2.8) 
below] 

(A fc NA 6 ) 1 ( (A b \b(iD±) 2 b\A b ) \ 1 ( (A b \ba ■ Gb\A b ) \ 

~ 1 + 2^f [ 2^ ) + 4^f { 2^ ) + ° [llm ^ (L2) 
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with D± = <9 M — v^v-D. This expression is not only model independent but also contains 
nonperturbative kinetic and chromomagnetic effects which are either absent or overlooked in 
the earlier quark-model calculations. 

Based on the OPE approach for the analysis of inclusive weak decays, predictions for the 
ratios of bottom hadron lifetimes have been made by several groups. The first correction to 
bottom hadron lifetimes is of order 1/mf and it is model independent [10]: 

T[in -l + 0(l/mg), 



r(B d ) 

r(Bs) - (1.00 ± 0.01) + Gil/ml), 



d 



t(B, 

^ = O.98 + 0(l/mg). (1.3) 

The I /ml corrections are small and essentially canceled out in the lifetime ratios. Nonspec- 
tator effects in inclusive decays due to the Pauli interference and ^-exchange contributions 
account for \/m\ corrections and they have two eminent features: First, the estimate of non- 
spectator effects is model dependent; the hadronic four-quark matrix elements are usually 
evaluated by assuming the factorization approximation for mesons and the quark model for 
baryons. Second, \/m\ corrections can be quite significant due to a phase-space enhancement 
by a factor of IQtc 2 . Predictions made in [11] for lifetime ratios of bottom hadrons are 

r(B d ) i - U + U - U5 l v 200MevJ ' r(B d ) ~ U * ' [LA) 

Experimentally [12], while the B~ and B d lifetimes are very close, it appears that the Ab 
lifetime is significantly shorter than the B meson one: 

= 1.06 ± 0.04 , = °- 79 ± °- 06 ( world average). (1.5) 

T{B d ) T{B d ) 

It should be mentioned that while the world average value for r(Ab) /r(B d ) is dominated by 
LEP experiments [12], the CDF experiment alone yields [13] 

4^4 = 0.87 ±0.11 (CDF). (1.6) 
T \ B d) 

It is thus important to fully settle down the experimental situation in the near future. Evi- 
dently, the conflict between experiment (1.5) and theoretical expectations from (1.3) or (1.4) 
is striking and intriguing. This has motivated several subsequent studies trying to understand 
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the enhancement of the A b decay rate [10,14-17]. For example, a model-independent analysis 
in [10] gives 



r(B d ) 
r(A 6 ) 



1 + 0.03^! + 0.004^2 - 0.70ei + 0.20e 2 , 

0.98 - 0.17ei + 0.20e 2 - (0.012 + 0.0215)r, (1.7) 



r(B d ) 

where £j, B^ B, r are the hadronic parameters to be introduced below in Sec. II. Note 
that while the ratio T(B~)/T(B d ) is predicted to be greater than unity in [11] [see (1.4)], 
it was argued in [10] that the unknown nonfactorizable contributions in (1.7) characterized 
by £j make it impossible to have reliable predictions on the magnitude of the lifetime ratio 
and even the sign of corrections. Since the measured ratio of t(B~) /r{B d ) is very close to 
unity, it follows from (1.7) that e\ « 0.3 e 2 [10]. Then it is clear that the data for the ratio 
r{Kb) /t{B(i) cannot be accommodated by the theoretical prediction (1.7) without invoking a 
too large value of r or B, which is expected to be order unity. It is reasonable to conclude 
that the 1/mf corrections in the heavy quark expansion do not suffice to describe the observed 
lifetime differences between A5 and B<±. 

In order to employ the OPE approach to compute inclusive weak decays of heavy hadrons, 
some sort of quark-hadron duality has to be assumed (for an extensive discussion of quark- 
hadron duality and its violation, see [18,19]). Consider the inclusive semileptonic decay. The 
OPE cannot be carried out on the physical cut in the complex v-q plane since T^, the time- 
ordered product of two currents, along the physical cut is dominated by physical intermediate 
hadron states which are nonperturbative in nature. To compute T^ v or the Wilson coefficients 
by perturbative QCD, the OPE has to be performed in the unphysical region far away from 
the physical cut. The question is then how to relate the operator product expansion for 
T^ u in the unphysical region to the physical quantities in the physical Minkowski space. 
Since the physically observable quantity is related to the imaginary part of T 7 ^, it can be 
reliably computed by deforming the contour of integration into the unphysical region [20,18], 
provided that the physical quantity involves certain integrals of T^ v in the physical region. 
This procedure is called "global duality" [18]. Global quark-hadron duality also means that 
the hadronic cross section is dual or matching to the OPE-based quark cross section. However, 
unlike the total cross section in e + e~ annihilation, there is a small portion of the contour near 
the physical cut where global duality can no longer be applied. As stressed in [18] , one must 
resort to local duality to justify the use of the OPE in this small region. Fortunately, the 



contribution is of order Aqcd / tuq and can be neglected for quantities smeared over an energy 
scale of order Aq C d- 

Global quark-hadron duality for inclusive semileptonic decays, namely the matching be- 
tween the hadronic and OPE-based expressions for decay widths or smeared spectra in semilep- 
tonic B and At, decays has been explicitly proved to the first two terms in l/mt, expansion 
and the first order in a s in the Shifman-Voloshin (SV) limit [21]. The hadronic decay rate is 
calculated by summing over all allowed exclusive decay channels. In the SV limit for B meson 
decays via b — > c transitions, the dominant hadronic final states are the D and D*. (At zero 
recoil, the quark-mixing-favored semileptonic decays of a B meson in the heavy quark limit 
can only produce a D or D* meson [21,22].) The exclusive decay rates or distributions for 
B — > (D + D*)£u depend on hadron masses, whereas the inclusive decay rates evaluated by the 
OPE depend on quark masses. Global duality is then proved by showing explicitly the equality 
of inclusive and exclusive decay rates. Note that this proof of global duality in QCD is valid 
only in the SV limit. Beyond this limit, it becomes difficult to sum over all allowed exclusive 
semileptonic decay channels and evaluate all of them. It was shown recently in [23] that a 
proof of quark-hadron global duality in the general kinematic region to order (Aq CD /7«b) 2 
can be achieved in the PQCD-based factorization approach, which is formulated in terms of 
meson-level kinematics rather than the quark-level one. It was demonstrated explicitly in [23] 
that the integrated quark-level spectrum equals to the hadron-level spectrum and that linear 
1/rrib corrections to the total decay rate are nontrivially canceled out, in agreement with the 
OPE expectation [7,8]. 1 

Unlike the semileptonic inclusive decays in which the use of the OPE is validated by 
deforming the contour away from the physical cut, it is pointed out in [18] that there is 
no external momentum q in inclusive nonleptonic decays which allows analytic continuation 
into the complex plane. Therefore, the OPE is a priori not justified in this case and local 
duality has to be invoked in order to apply the OPE directly in the physical region. It is 
obvious that local quark-hadron duality is less firm and secure than global duality, although 
its validity has been proved to the first two terms in 1/mg expansion and the first order in 



1 The absence of linear l/m^ corrections to decay widths is trivial in the SV limit since the inclusive 
decay rates depend on AM = tub — mp rather than mg, and AM = 5m + 0(1/ ml) with 5m = 
m b - m c [21]. 
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a s in the SV limit under the factorization hypothesis [24]. It should be stressed that quark- 
hadron duality is exact in the heavy quark limit, but its systematical 1/mg expansion is still 
lacking. It is very likely that 1 / ijiq corrections to quark-hadron duality behave differently for 
inclusive semileptonic and nonleptonic decays. Motivated by the conflict between theory and 
experiment for the lifetime ratio r(A(,)/r(Sd), it was suggested in [16] that the assumption 
of local duality is not correct for nonleptonic inclusive width and that the presence of linear 
l/rrib corrections is strongly indicated by the data. Moreover, the l/m& corrections are well 
described by the simple ansatz that the heavy quark mass ijiq is replaced by the decaying 
hadron mass in the m b q factor in front of all nonleptonic widths. It is easily seen that the 
factor (m B /mA t ) 5 = 0.73 is very close to the observed value of T(A b )/r(B d ) . Under this 
ansatz, a much better description of lifetimes of bottom and charmed hadrons was shown in 
[16]. Irrespective of the lifetime ratio problem, there is another important reason why this 
ansatz is welcome. The absolute decay rate of the B meson predicted in the OPE approach 
is at least 20% smaller than the experimental value (see Sec. Ill below). We shall show in 
Sec. Ill that the discrepancy between theory and experiment is greatly improved when the 
nonleptonic width scales with m%. 

In the aforementioned factorization approach of [23] , the nonleptonic width of bottom 
hadrons scales with m 5 Hb . Local duality means that a replacement of meson- level kinematics 
by quark kinematics, for example, rriH b = rrib(l + A# 6 /m& + •••),••• etc., will turn T^l into 
r^L E , the OPE-based decay rate. Consequently, the relation between violation of local duality 
and the above-mentioned ansatz will become natural in the factorization approach. 

In the present paper we will study nonspectator effects in inclusive nonleptonic and semilep- 
tonic decays and analyze the lifetime pattern of heavy hadrons. In particular, we focus on 
the lifetimes of heavy baryons and study the implications of broken local duality. We will 
demonstrate that the lifetime hierarchy of bottom baryons is dramatically modified when the 
quark mass is replaced by the hadron mass in nonleptonic widths. The layout of this paper is 
organized as follows. In Sec. II we give general heavy quark expansion expressions for inclusive 
nonleptonic and semileptonic widths and pay attention to the evaluation of baryon four-quark 
matrix elements and the nonperturbative parameter A2 for baryons. We then study bottom- 
hadron lifetimes in Sec. Ill and apply the ansatz mentioned above. In Sec. IV we examine the 
applicability of the same prescription to charmed baryon decays. Discussions and conclusions 
are given in Sec. V. 
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II. FRAMEWORK 



In this section we write down the general expressions for the inclusive decay widths of heavy 
hadrons and evaluate the relevant hadronic matrix elements. It is known that the inclusive 
decay rate is governed by the imaginary part of an effective nonlocal forward transition oper- 
ator T. When the energy released in the decay is large enough, the nonlocal effective action 
can be recast as an infinite series of local operators with coefficients containing inverse powers 
of the heavy quark mass rag. Under this heavy quark expansion, the inclusive nonleptonic 
decay rate of a heavy hadron Hq containing a heavy quark Q is given by [7,8] 

(I o (x,0,0) + I o (x,x,0)){H Q \QQ\H Q ) 
1 



4 2cic 2 



2 (h(x, 0, 0) + h(x, x, 0))(H Q \Qa ■ GQ\H Q ) 

(I 2 (x, 0, 0) + I 2 (x, x, 0))(H Q \Qa ■ GQ\H Q ) 



m Q N c 

+-J—(H Q \C nspec \H Q ) + 0{l/m%), (2.1) 
2m HQ 

where a ■ G = o^G^ , Ci, c 2 are Wilson coefficient functions, N c = 3 is the number of color, 
the factor £ takes care of the relevant CKM matrix elements, for example, £ = | Vrf,V^j 2 for 
quark- mixing-favored bottom decay, Jo, h and I 2 are phase-space factors: 

I (x, 0, 0) = (1 - x 2 )(l -8x + x 2 ) - 12x 2 lnx, 
h(x, 0, 0) = i(2 - x^)I (x, 0, 0) = (1 - x)\ 

/ 2 (x,0,0) = (l-xf (2.2) 

for b — > cud (x = m 2 /ml) or c — > sud (x = m 2 /m 2 ) transition and 

I (x, x, 0) = v(l - Ux - 2x 2 - 12x 3 ) + 2Ax 2 {\ - x 2 ) In Ltl, 

1 — v 

h(x,x,G) = ^(2- x-^-)I (x,x, 0), 

I 2 (x, x, 0) = v(l + ^ + 3x 2 ) - 3x(l - 2x 2 ) In (2.3) 

for b — > ccs transition with t> = \/l — 4x. 

The dimension-six four-quark operators £ nS pcc in (2.1) describe nonspectator effects in 
inclusive decays of heavy hadrons and are given by [3-5] 



£nspec = £ (1 - X) 2 {{ C \ + C^) (QQ) {Ml) + 2 Ci C 2 (QQi) (QlQ) } 



2tt 

9. f 

X, 



1 1 " ln dc\{l - xf[{l + X -)(QQ)(M2) - (l + 2x)Q a (l - 75 )gM(l + 7 5 )Q1 



67T 



+ (2 Cl c 2 + N c c\){l - x) 2 [(l + |)(Qg 2 )(g 2 Q) - (1 + 2x)Q(l - 75)9292(1 + 7 5 )Q] j 
- ejc? v/T^4i [(1 - x)(QQ)(«,(b) - (1 + 2*)Q a (l - 7s)«?^(l + 7 5 )Q Q ] 

+ (2 Cl c 2 + 7V c ^)v / T 3 4^ [(1 - x)(Qq 3 )(q 3 Q) - (1 + 2x)Q(l - 75)^3(1 + 7s)<5] 
67 J<|(1 - ^) 2 [(1 + ?)(QQ)(«i(fe) - (1 + 2x)Q a (l - 75)^^(1 + 75)Q1 



67T 

+ (2 Cl c 2 + N c c\){\ - x) 2 [(l + ^)(0g 3 )feQ) - (1 + 2x)Q(l - 75)9353(1 + 7s)Q] j, (2.4) 

where (g'g) = g'7 M (l — 75)9, and a, j3 are color indices. Note that for charm decay, Q = 
c, qi — u, q-i — d and q 3 = s and for bottom decay, Q = b, q 1 = u, q 2 = d, q 3 = s. 
The last term in (2.4) is due to the constructive interference of the s quark and hence it 
occurs only in charmed baryon decays. The third term in (2.4) exists only in bottom decays 
with cc intermediate states. For inclusive nonleptonic decays of heavy mesons, the first term 
in (2.4) corresponds to a Pauli interference and the second and third terms to W^-exchange 
contributions. For heavy baryon decays, the first term is a W-exchange contribution and the 



rest are interference terms. The phase-space suppression factors e.g. (1 — x) 2 , y/l — Ax, ■ ■ ■ 
etc. in (2.4) are derived in [25,10]. 

Several remarks are in order, (i) There is no linear I/tuq corrections to the inclusive decay 
rate due to the lack of gauge-invariant dimension-four operators [20,7], a consequence known 
as Luke's theorem [26]. Nonperturbative corrections start at order 1/toq. (ii) It is clear 
from Eqs. (2.1) and (2.4) that there is a two-body phase-space enhancement factor of IQtt 2 
for nonspectator effects relative to the three-body phase space for heavy quark decay. This 
implies that nonspectator effects, being of order 1/mg, are comparable to and even exceed the 
1/rriQ terms, (iii) For charmed meson decay, the 1/N C correction to T NL characterized by the 
term (2ciC2/N c ){H c \cc\H c ) is found to be compensated by the nonperturbative gluonic effect 
[i.e. the term proportional to I 2 (x, 0, 0)]. This cancellation is small for B meson decay due to 
the smallness of \/m\. This indicates that the rule of discarding 1/N C terms [27] is operative 
in charm decays but not so for the B meson case, (iv) Thus far the Wilson coefficients and 
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four-quark operators in Eq. (2.4) are renormalized at the heavy quark mass scale. Sometimes 
the so-called hybrid renormalization [5,28] is performed to evolve the four-quark operators 
(not the Wilson coefficients !) from mg down to a low energy scale, say, a typical hadronic 
scale yUhad- The underlying reason is that the factorizable approximation for meson matrix 
elements and the quark model for baryon matrix elements are believed to be more reliable at 
the scale /Zhad- The evolution from mq down to /ihad will in general introduce new structures 
such as penguin operators. However, in the present paper we will follow [10] to employ (2.1) 
and (2.4) as our starting point for describing inclusive weak decays since it is equivalent to 
first evaluating the four-quark matrix elements renormalized at the uiq scale and then relating 
them to the hadronic matrix elements renormalized at /ihad through the renormalization group 
equation, provided that the effect of penguin operators is neglected. 

For inclusive semileptonic decays, apart from the heavy quark decay contribution there is 
an additional nonspectator effect in charmed-baryon semileptonic decay originating from the 
Pauli interference of the s quark [29]. It is now ready to deduce the inclusive semileptonic 
widths from (2.1) and the last term in (2.4) by putting c\ — 1, c 2 = and N c = 1: 



G^m^ 2 r)(x,x e ,0) 



1 



T Sh (H Q ) = -L-^-\Vckm\ 2 J o(^ °> 0){Hq\QQ\H q ) - — h(x, 0, 0)(H Q \Qv ■ GQ\H Q ) 



1927T 3 ' ™' 2m Hq 



x 



. XcsV- (l-x) 2 (l + -)(cs)(sc)-(l + 2a;)c(l-75)ss(l + 75)c, (2.5) 

07T 2rriH c 1 2 J 

where rj(x, X£,0) with x e = {m^/mQ) 2 is the QCD radiative correction to the semileptonic 
decay rate. Its general analytic expression is given in [30]. The special case r](x, 0, 0) is given 
in [31] and it can be approximated numerically by [32]: 



^,0,0)- 1-^ 



:- 2 -^)(i-v^) 2 + ^ 



(2.6) 



With x = and the replacement a s — > |a, (2.6) is reduced to the well-known QED correction 
to the muon decay. The second term in Eq. (2.5) occurs only in the semileptonic decay of E c 
and Q c baryons. 

We next turn to the 2-body matrix elements (Hq\QQ\Hq). The use of the equation of 
motion 

QQ = Qj)Q + -^QitD^Q + -^Qa ■ GQ + 0(\/m%), (2.7) 

ZTUq 41TLq 

with Dl = - v^v-D, leads to 



(H Q \QQ\H Q ) _ K H G H 
2m Hg 2m Q H' 

with 

K H = --^(H Q \Q(iD ± ) 2 Q\H Q ) = -X u 
Zm Hq 

G H = ^—(HqIqI* ■ G Q\ R q) = ( 2 - 9 ) 
The mass of the heavy hadron Hq is then of the form 

m HQ = m Q + A HQ - — -—, (2.10) 

where the three nonperturbative HQET parameters A# Q , Ai, A2 are independent of the heavy 
quark mass and in general Ah q is different for different heavy hadrons. Since a ■ G ~ Sq ■ B 
and since the chromomagnetic field is produced by the light cloud inside the heavy hadron, it 

— * — * — * — * 

is clear that a ■ G is proportional to Sq ■ Si, where Sq (Se) is the spin operator of the heavy 
quark (light cloud). More precisely, 

d H = -{H Q \4S Q .S t \H Q ) 

= -2[S tot (Stot + 1) - Sq(S q + 1) - S e (S e + 1)]. (2.11) 

Therefore, d H — 3 for B, D mesons, d H — — 1 for B*, D* mesons, d H = for the antitriplet 
baryon Tq, du = 4 for the spin-| sextet baryon Sq and dn = —2 for the spin-| sextet baryon 
S Q . It follows from (2.10) that 

x meson _ 1 / 2 m 2 x _ / 0.12 GeV 2 for B meson; 

An — —[flip* — fflpl — S r> 

2 4 V p PJ \ 0.14 GeV for D meson, 

\r^ = \{ml* Q -ml Q ). (2.12) 

The values of A^ 011 will be fixed later. As for the kinetic energy parameter Ai we use [33] 

A mcson _ A baryon = ± Q ^ ^2 _ 

This leads to 



m b -m c = ((m B ) - (m D )) 1 



2{m B ){m D ) J 

(3.40 ± 0.03) GeV, (2.14) 
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where (mp) = \(m P + 3m P ») denotes the spin-averaged meson mass. 

We will follow [10] to parametrize the hadronic matrix elements in a model-independent 
way. For meson matrix elements of four-quark operators, we follow [10] to define the param- 
eters Bi and ef 

(B q \(bq)(qb)\B q ) = f 2 B m 2 B B u 

(B q \b(l - 75)99(1 + l 5 )b\B q ) = f 2 Bq m 2 Bq B 2 , 

(B q \(bt a q)(qt%)\B q ) = f B m 2 Bq e u 

(B q \bt a (l - lb )qqt a (l + lb )b\B q ) = f 2 B m 2 Bq e 2 , (2.15) 

where (q't a q) = g't a 7 M (l — 75)^ and t a = X a /2. Under the factorization approximation, Bi 
and are given by Bi = 1 and = 0, but they will be treated as free parameters here. As a 
consequence of (2.15), we obtain 

(B g \(bb)(qq)\B g ) = f B m\ q + 2 £l ) , 

(^1^(1 - ^^(l + 75)fc Q |S ? ) = flm\ {^-B 2 + 2£ 2 ) . (2.16) 

As for the baryon matrix elements of four-quark operators we have to rely on the quark 
model. We first consider the MIT bag model [34] and define three four-quark overlap integrals: 

a q = J d 3 r[u 2 q (r)u 2 Q (r) + v 2 (r)v 2 Q (r)], 
b q = /rf 3 r[ M J(r)^(r)+^(r)^(r)], 

Cq = J d 3 ru q (r)v q (r)u Q (r)v Q (r), (2.17) 

which are expressed in terms of the large and small components u(r) and v(r), respectively, of 
the quark wavefunction. For the antitriplet heavy baryon Tq or the sextet heavy baryon Qq 
(recall that only the 0° and il^ of the sextet baryons decay weakly), the four baryon matrix 
elements 

(T Q I (Qq) (qQ) \T Q ) , (Tq \ (QQ) (qq) \T Q ) , 

(T Q \Q(l - 75)99(1 + 1*)Q\T Q ), (T Q \Q a (l - 75 )9V(1 + ls)Q a \T Q ) 
are not all independent. First of all, we have 

(T Q \(Qq)(qQ)\T Q ) = -(a q + b q )(2m TQ ), 

(Q Q \(Qs)(sQ)\n Q ) = ~(l8a s + 2b s + 32c s )(2m nQ ), (2.18) 
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(see e.g., Ref. [35] for the technical detail of the bag model evaluation), where we have taken 
into account the fact that there are two valence s quarks in the wavefunction of the Qq. 
Second, since the color wavefunction for a baryon is totally antisymmetric, the matrix element 
of (QQ)(qq) is the same as that of (Qq)(qQ) except for a sign difference. Thus we follow [10] 
to define a parameter B 

(T Q \(QQ)(qq)\T Q ) = -B{T Q \(Qq)(qQ)\T Q ), 

(Q Q \(QQ)(ss)\Q Q ) = -B{Q Q \(Qs)(sQ)\Q Q ), (2.19) 
so that B — 1 in the valence-quark approximation. Third, it is straightforward to show that 
(W q 7m7 5 Q / W(1 - 7 5 )<f |T Q > = 0, 

(^|Q q 7m75Q / W(1 " 7 5 )<f |^q> = 4 (a - ^ (2mn Q ). (2.20) 

The first relation in (2.20) is actually a model-independent consequence of heavy-quark spin 
symmetry [10]. Since 

Q^TsQVra - 7 5 )g a = -0(1 - 75)55(1 + 75 )Q - \(Qq)(qQ), (2.21) 
it follows from (2.20) that 

(T Q \Q a (l - 7s)5V(l + K)Q a \T Q ) = -B{T Q \Q{1 - 75)55(1 + 1 5 )Q\T Q ) 

= -- J B(a 9 + 6 g )(2m TQ ), 

(Q Q \Q a (l - 75 )sV(l + l5 )Q a \tt Q ) = -B(Q Q \Q(1 - 75 )ss(l + l5 )Q\Q Q ) 

5 16 

= B(a s --b s --c s )(2m nQ ). (2.22) 

In the nonrelativistic quark model (NQM), baryon matrix elements of four-quark operators 
are the same as that of (2.18) and (2.22) except for the replacement: 

a q - l^(0)| 2 = / d 3 ru 2 q (r)u 2 Q (r), b q ^0, c q ^ 0. (2.23) 

In general, the strength of destructive Pauli interference and ^-exchange is governed by 
a q + b q in the bag model and |^(0)| 2 in the NQM. However, it is well known in hyperon decay 
that the bag model calculation of a q + b q gives a much smaller value than the nonrelativistic 
estimate of |^(0)| 2 : a u + b u ~ 3 x 10- 3 GeV 3 , while |^(0)| 2 ~ 10- 2 GeV 3 . We shall see later 
that this also occurs in bottom baryon decay. As pointed out in [36], naively one may be 
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tempted to conclude that the relativistic models are presumably more reliable. For example, 
the lower component of the wavefunction is needed to reduce the NQM prediction g& — § to 
the experimental value of 1.25. However, the difference between a u + b u and |-0(O) | 2 is not 
simply attributed to relativistic corrections; it arises essentially from the distinction in the 
spatial scale of the wavefunction especially at the origin. As a consequence, both models give 
a quite different quantitative description for processes sensitive to |-^(0)| 2 . It has been long 
advocated in [37] that a small value of |-0(O) | 2 should be discarded since a realistic potential 
that fits to the orbital-excitation spectrum yields (S(f 1 — f 2 )) ~ 10~ 2 GeV 3 . Empirically, it 
also appears that the NQM works better for charmed baryon decays [4,36]. 

In the following we will consider the NQM estimate of baryon matrix elements. Consider 
|^ 6 (0)| 2 as an example. A straightforward calculation of hyperfine splitting between S fe and 
A b yields [38] 

m^ b -m Ab = — « S N- rH< 6 (0)| 2 , (2.24) 

9 m^m^ 

where the equality |Vv(0)| 2 = |Vv(0)| 2 has b een assumed. The uncertainties in Eq. (2.24) 
associated with a s {m h ) and the constituent quark mass m q can be reduced by introducing the 
.B-meson wavefunction at the origin squared |^(0)| 2 = -^f B m B which is related to the B* 
and B mass difference by m B * — rn B = Y7ro; s (mft)|^(0)| 2 /(m b m g ). Hence, 

l< 6 (0)| 2 = |^_ (0)|2 _ (2 . 25) 

y m b — m q m B * — m B 

Another method is proposed by Rosner [15] to consider the hyperfine splittings of and B 
separately so that 



K 6 (o)i 2 = \^m\ 2 = \ m I b r* '< (0)|2 - (2 - 26) 

O II l B * lli B 



This method is supposed to be most reliable as |^6 g (0)| 2 thus determined does not depend on 
a s and m q directly. Numerically, we find that (2.25) and (2.26) both give very similar results. 
Defining the wavefunction ratio 

2 

(2.27) 



< 6 (0) 



<(0) 

the baryon matrix elements in (2.18) and (2.22) can be recast to 
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(T b \(bb)(qq)\T b ) = - B (T b \(bq) (qb)\T b ) = ±f B m Bq rB(2m Tb ), 
(T b \b(l - 75)55(1 + l5 )b\T b ) = ^f Bq m Bq r(2m Tb ), 
(T b \b a (l- l5 )qPf(l + l5 )b a \T b ) = -±f B m Bq rB(2m Tb ), 
(Q b \(bb)(ss)\Q b ) = -B(Q b \(bs)(sb)\Q b ) = ±f B m Bq rB(2mn b ), 
(Q b \b(l - 75)^(1 + l5 )b\n b ) = -^f 2 Bq m Bq r(2mn b ), 

(n b \b a (l - 75)^(1 + l5 )b a \il b ) = ^f 2 Bq m Bq rB(2mn b ), (2.28) 

where f Bq is the decay constant of the meson B q . 

To estimate \ipbq(0) | 2 and the parameter r in the NQM, we find from (2.26) 2 

_ 4 m n - m Sh _ 4 m B; - m B > _ 4 m n; - to^ 

3 to#. — rng 3 m B * — m B 3 m B * — m B 

and likewise for r^ c , r= c and r^ c , where E' bc denote spin-| sextets. Heavy baryon masses 
have been studied in [39] in 1/mg and l/N c expansions within the HQET framework. The 
chromomagnetic mass splittings for charmed baryons are given by [39] 

m E . - m Ec = 65.7 ± 2.3 MeV, to h * - m E , c = 63.2 ± 2.6 MeV, 

msi, - m Qc = 60.6 ± 5.7 MeV, (2.30) 

where precise measurements of S* and H* have been reported by CLEO [40]. It is evident 
that the heavy-quark spin-violating mass relation [39] 

(m E * - to Sc ) + (toq* - toqJ = 2(to h » - to H /J (2.31) 

is very accurate. It follows that 




(2.32) 



for m b = 5 GeV, m c = 1.6 GeV (see below). This mass splitting is substantially smaller than 
the preliminary result to s * — to S(; = (56 ± 16)MeV reported by the DELPHI Collaboration 



2 Our result for r\ b is the same as [15] but different from [10] in which r\ h is given by |(m s * — 

m l b )/( m B* -™>%). 
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[41]. Since the measured mass difference of S* and S c is around 66 MeV [cf. (2.30)], a large 
hyperfine splitting of order 55 MeV for the baryon is very unlikely. Likewise, 

m E * -m s > = 20.2 MeV, m Q * - ma = 19.4 MeV. (2.33) 

b ' 'b b 

Because Am B = m B * - m B = 45.7 ± 0.4 MeV and Am D = m D * - m D = 143 MeV [42] [note 
that ArriB and Amu obey the same scaling relation as (2.32)], we find 



r Ac = r Ab = 0.61, r Hc r Sb = 0.59, r Qc m b = 0.53, (2.34) 



and 



|^(0)| 2 = 0.87 x 10~ 2 GeV 3 , |$*(0)| 2 = 0.84 x lO^GeV 3 , 

|^(0) | 2 = 0.81 x 10~ 2 GeV 3 , (2.35) 

for f Bq = 180 MeV [43]. An estimate in the QCD sum rule analysis yields r ~ 0.1 — 0.3 
[17]. Therefore, the NQM estimate of 1^(0) | 2 is indeed larger than the analogous bag model 
quantity: a q + b q ~ 3 x 10~ 3 GeV 3 . However, for the charmed baryon we obtain |V^ C (0)| 2 = 
3.8 x 10 _3 GeV 3 for fp = 200 MeV [43], which is smaller than those in bottom or hyperon 
decay. It seems that the smallness of IV^WI 2 * s ascribed to the assumption that the D meson 
wavefunction at the origin squared |^(0) | 2 is given by -^f D m D . We will come back to this 
point in Sec. IV. By comparing (2.28) with (2.18) we see that r is of order 0.20 in the bag 
model. 

Finally we are ready to estimate the HQET parameter js^ 1117011 [see Eq. (2.12)]. Using the 
baryon masses [39] 

m Sc = 2452.9 MeV, m^ c = 2580.8 MeV, m Qc = 2699.9 MeV, 

m T, b = 5824.2 MeV, m s > = 5950.9 MeV, m^ b = 6068.7 MeV, (2.36) 

and (2.30)-(2.33) we find 

0.055 GeV 2 for charmed baryons; 



^baryon _ ^ 



0.041 GeV 2 for S 6 ; 
0.040 GeV 2 for E' b ; 



U.U4U KjC \ 

k 0.039 GeV 2 for Q b . 



(2.37) 



It is interesting to note that the large- N c relation [39] 

A mcson _ jy^baryon ^ ^ 

is fairly satisfied especially for bottom hadrons. 
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III. LIFETIMES OF BOTTOM HADRONS 



Using the formulism described in the last section, semileptonic and nonleptonic widths 
are calculated in this section. We shall first try to fix the heavy quark pole mass from the 
measured inclusive semileptonic decay rate. The semileptonic width of the B meson given by 
(2.5) 

rsl( B^ M )= ^|^,0,0) {l„(»,(^0)M-^(»,0,0)<^^} (3.1) 
1927r d [ 2m B m b 2m B J 

has the salient feature that empirically Tsl(B) is very insensitive to the choice of m& as long 
as m h — m c , which is free of renormalon ambiguity is fixed according to Eq. (2.14). Hence, 
we may use the measured T SL (D) to fix m c to be 1.6 GeV (see Sec. IV below) which in turn 
implies mj = 5 GeV, in excellent agreement with the pole mass determined from lattice QCD: 
mj = 5.0 ± 0.2 GeV [44]. Since x = (m c /mb) 2 = 0.1024, the phase-space factors U in (2.2) 
and (2.3) read 

I (x, 0, 0) = 0.476, l (x,x,0) = 0.147, h(x, 0, 0) = 0.649, 

h(x,x,0) = 0.328, I 2 (x, 0,0) = 0.723, I 2 (x, x, 0) = 0.220. (3.2) 

From (3.1) we obtain 

T(B -> Xev) = AAA x 10" 14 GeV, 

T SL (B) = 2.24 T(B -> Xev) = 9.95 x 10- 14 GeV, (3.3) 
for \ Vcb\ = 0.039, where uses of Eqs. (2.8), (2.9), (2.12) and (2.13) have been made, for example, 

6A meso„ = Q 72 Ge y2_ (3 ^ 



'\B\ba-Gb\B) ,, XJ1K , (1 ^ n -,, r ,. A - 2 



2m B 

Since the phase space for the r semileptonic decay mode relative to that of the e mode is 
0.24 : 1, this accounts for the factor 2.24 in (3.3). The result (3.3) agrees very well with 
experiment [42] 

T(B-/B° admixture -> Xev) = (4.31 ± 0.17) x 10~ 14 GeV. (3.5) 

Likewise, we find for bottom-baryon semileptonic decays 

T(A 6 -> Xev) = T(E b -> Xez/) = 4.59 x 10~ 14 GeV, 

T(Q b -»• XeP) = 4.53 x 10~ 14 GeV, (3.6) 
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and hence 

r SL (A b ) = r SL (S 6 ) = 2.24 T(A 6 -> Xev) = 1.027 x 10- 13 GeV, 

T SL (n b ) = 2.2AT(n b -> ATeP) = 1.014 x HT 13 GeV. (3.7) 

Note that the tiny difference between r SL (A b ) and r gL (fif,) arises from the fact that the 
chromomagnetic operator contributes to the matrix element of Q b but not to A b (or E b ) as 
the light degrees of freedom in the latter are spinless; that is, 



{Ablla ' GblAb) = 0, WhT ' GbW = 8Ar y ° n = 0.31 GeV 2 . (3.8) 
2m Afc 2m Qb 

To compute the nonleptonic decay rate we apply the Wilson coefficient functions 

d(^) = 1.14, c 2 (/i) = -0.31, (3.9) 

which are evaluated at fi — 4.4 GeV to the leading logarithmic approximation (see Table XIII 
of [45]). From Eq. (2.1) the nonleptonic widths of bottom baryons arising from b quark decay 
are found to be 

T dec (B) = 2.216 x 10~ 13 GeV, T dec (tt b ) = 2.217 x 10~ 13 GeV, 

r dec (A 6 ) = T dcc (S 6 ) = 2.220 x 10- 13 GeV. (3.10) 

We see that the b quark decay contribution r dcc is very similar for bottom hadrons even 
though the chromomagnetic mass splitting is different among them. Therefore, to C(l/m 3 ) 
we obtain 

We next turn to the nonspectator effects of order 1/m 3 . The Pauli interference in inclusive 
nonleptonic B~ decay and the VF-exchange contribution to B d can be evaluated from the first 
and second terms in (2.4): 



1 



r*™(B d ) = -r Vnspcc \ (1 - x)\l + -x) (-of + 2 Cl c 2 + N c c I 2 )B l + 2c{e x 



1 



3 



-(1 - xfd + 2x) [{^c\ + 2 Cl c 2 + N c c\)B 2 + 2c\e 2 ] |, 
r?(B~) = r ?7n S pec(l - xf [{c\ + c 2 2 )(B 1 + 6 £l ) + Qc^B^ , (3.12) 



with [10] 
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r = S4l^! 2 ' %spcc = 167T 2 ^, (3.13) 

where we have applied Eqs. (2.15), (2.16) and neglected Cabibbo-suppressed ^-exchange 
contribution to B d . As stressed in [10], the coefficients of Bi in (3.12) are one to two orders of 
magnitude smaller than that of E{. Therefore, the contributions of Bi can be safely neglected 
at least in T ann (B d ). Numerically, 

T ann (B d ) = (-0.49l£i + 0.563e 2 ) x l(T 13 GeV, 

T^(B-) = (-0.130Bi + 1.505ei) x 10~ 13 GeV. (3.14) 

Beyond the factorization approximation, £j may receive nonfactorizable contributions. A QCD 
sum rule estimate gives S\ ~ —0.15 and e 2 ~ [25]. This implies a constructive W-exchange 
to B d and a destructive Pauli interference to B~ . 

As for the nonspectator effects in nonleptonic decays of bottom baryons we obtain from 
(2.4) that 



T ann (A 6 ) = r 77nspcc r(l - x) 2 (B{c\ + 4) - 2 Cl c 2 ), 



r: it (A 6 ) = -ir r ?nspec r 



'l-a;) 2 (l + a;) + 



V c d 2 



V, 



ud 



y/l-Ax 



(Bc{-2c 1 c 2 -N c c 2 2 ) ) 



r ann (H°) = r ann (A 6 ), r hlt (s°) = r™\A b ), 

r^(S 6 -) = -Ir 7? nspcc r[(l - xf{\ + x) + VT^4x~] (Bc\ - 2 Cl c 2 - N c c 2 2 ), 
T hlt (fi 6 -) = -ir r ?nspcc rv / T^4i(5 - %x)(Bc{ - 2 Cl c 2 - N c cj), (3.15) 



4 

.1] 

6 

where use has been made of (2.28). Note that there is no W^-exchange contribution to the 



2^ and Qb and that there are two Cabibbo-allowed Pauli interference terms in decay, and 
one Cabibbo-allowed as well as one Cabibbo-suppressed interferences in A& decay. It is easily 
seen that under the valence-quark approximation i.e. B — 1, the W^-exchange contribution 
pann j g p r0 p 0r tional to c_ = (ci — c 2 )/2 as the four-quark operator + = {(jiq 2 ) {(jsQA) + 
(Q1Q4) (9392) is symmetric in color indices whereas the color wavefunction for a baryon is totally 
antisymmetric. Writing 

r NL = r dcc + r ann + r to , (3.16) 

the numerical results for nonleptonic inclusive decay rates are 
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2.220 + (0.042 + 0.0585)r] x 10~ 13 GeV, 
2.220 + (0.043 + 0.0665)7-] x 10~ 13 GeV, 
2.220- (0.037 + 0.1145)^ x 10~ 13 GeV, 
2.217- (0.043 + 0.1335)rj x 10~ 13 GeV, 



(3.17) 



where for later convenience we have normalized the parameter r in (3.17) to r\ b [see (2.34)]; 
that is, we have taken into account SU(3) breaking effect for r. Note that Bi in (3.14) and 
B, r in (3.17) are all renormalized at fx = 4.4 GeV. 

Before proceeding, it is worth emphasizing the difference between the W^-exchange con- 
tributions in the inclusive nonleptonic decays of the B meson and the bottom baryon. It is 
conventionally argued that IT-exchange in heavy meson decay is suppressed by helicity and 
color mismatch. For example, ^-exchange in B decay is helicity suppressed by a factor of 
167T 2 (/s/ms) 2 relative to the heavy quark decay amplitude. 3 By contrast, W^-exchange in 
baryon decay is neither helicity nor color suppressed. The diquark Qq system in the heavy 
baryon can have a spin configuration and the decay of a spin (not spin 1 !) state into two 
quarks is not subject to helicity suppression. 

Since B is of order unity and r ~ 0.60, it is evident from (3.17) and (3.10) that the bottom 
baryon lifetimes follow the pattern (see also Table I below) 



This pattern originates from the fact that while S^~, fi& all receive contributions from 

destructive Pauli interference, only A b and S° have ^-exchange and that r mt is most large 
in fl b due to the presence of two valence s quarks in its quark content. We shall see shortly 
that this lifetime pattern is dramatically modified when the b quark mass is replaced by the 
bottom baryon mass in nonleptonic widths. 



3 It had been claimed that soft gluon emission from the initial quark line or soft gluon content 
in the initial wavefunction can vitiate both helicity and color suppression [46]. The net effect is 
that the factor fs/f^B is effectively replaced by /g/m,, where m q is the constituent quark mass 
of the antiquark in the B meson [47]. As a consequence, contributions of TT-exchange will exhibit 
powerlike (mB/m q ) 2 enhancement and this renders the treatment of the heavy quark expansion 
for IT-exchange invalid. This issue was resolved by Bigi and Uraltsev [47] who showed that such 
powerlike enhancement does not arise for fully inclusive transitions and the soft gluon effect merely 
amounts to renormalizing the coefficients of 4-quark operators. 



r(n>-) ~ r(S 6 -) > r(A°) ~ r(3g). 



(3.18) 
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It follows from (3.3), (3.7), (3.10), (3.14) and (3.17) that 

IiM = 0.99 - 0.15ei + 0.17^ - (0.013 + 0.018S)r, (3.19) 

T \ B d) 

which is a model-independent result. This is consistent with the result (1.7) obtained in 
[10] with £j, B, r renormalized at fx — 4.85 GeV and with f B = 200 MeV. As stated in 
the Introduction, S\ and e<i obey the constraint S\ ~ 0.3 £2, then it is quite difficult, if not 
impossible, to accommodate the experimental value (1.5) for r(A 6 )/r(S d ) without invoking too 

large value of r and/or B. We will argue below that the contribution of — 0.15eH 0.018i?r 

in (3.19) is at most of order 6%. We hasten to remark that the current CDF result (1.6) for 
the lifetime ratio is consistent with the theoretical prediction (see, however, a comment after 
Eq. (3.22)). 

Irrespective of the above-mentioned lifetime ratio problem, there exists another serious 
difficulty, namely the predicted absolute decay width of the B or A b hadron based on the heavy 
quark expansion [see (3.3), (3.7), (3.10) and (3.17)] is too small compared to the experimental 
values [12]: 

F(B d ) = (4.2461°;°^) * 10" 13 GeV, r(B d ) = (1.55 ± 0.04) ps, 
T(B-) = (3.965l^|) x 10~ 13 GeV, t(B~) = (1.66 ± 0.04) ps, 
r(A 6 ) = (5.351^|) x 10~ 13 GeV, r(A 6 ) = (1.23 ± 0.09) ps. (3.20) 

Obviously, even if the destructive contribution r mt (i?~) is not taken into account, the result 
T dcc {B) + r SL (fi) = 3.211 x 10~ 13 GeV is too small by about 20% to account for the observed 
decay rate of B~ . 4 To compute the decay widths of bottom baryons, we have to specify the 



4 The problem with the absolute total decay width T(B) of the B meson is intimately related to 
the problem with the 5-meson semileptonic branching ratio -Bsl- The theoretical prediction for -Bsl 
is in general above 12.5% [48], while experimentally -Bsl = (10.23 ± 0.39)% [49]. In our case we 
obtain Bsl ^ 13.8%. Several scenarios have been put forward in the past to resolve the discrepancy 
between theory and experiment for Bsl or T(B). Here we mention two of the possibilities, (i) Since 
the theoretical results depend on the scale li to renormalize a s (/f) and the Wilson coefficients c\^{^), 
one may choose a low renormalization scale, n/rrn, ~ 0.3 — 0.5, to accommodate the data [10]. Local 
duality holds in this scenario, (ii) Next-to-leading order QCD radiative corrections to nonleptonic 
decay will increase the rate for b — > ccs substantially and decrease -Bsl [50,51]. Using the result 
of [50], we find that the QCD effect will bring Bsl down by 1% and hence Bsl ^ 12.7%. It was 
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values of B and r. Since B = 1 in the valence-quark approximation and since the wavefunction 
squared ratio r is evaluated using the quark model, it is reasonable to assume that the NQM 
and the valence-quark approximation are most reliable when the baryon matrix elements 
are evaluated at a typical hadronic scale //had- As shown in [10], the parameters B and r 
renormalized at two different scales are related via the renormalization group equation to be 



B(n)r(ii) = J B(/i had )r(/i had ) 



(3.21) 



with 



K 



had J 



3JV c /2/3 



a s (n) 



a s ((J> 



had J 



(3.22) 



and p — — |n/. Choosing a s (/ihad) = 0.5 and /i = 4.4 GeV, we obtain B(/jl) = 

0.59B(fi had ) ~ 0.59 and r(fx) ~ 1.7r(fi had ). Using r(fi had ) = 0.61 [see (2.34)], the calculated 
decay rates of bottom baryons are summarized in Table I. It is evident that the predicted A b 
lifetime is too large by 8 standard deviations. Note that while the CDF measurement (1.6) 
for the lifetime ratio r(Ab)/r(5d) can be easily accommodated in theory, it is still difficult to 
explain the absolute lifetime r(A 6 ) = (1.32 ± 0.15 ± 0.07)ps measured by CDF [13]. 



Table I. Various contributions to the decay rates (in units of 10 GeV) of bottom baryons 





pdcc 


pann 


pint 


r SL 


ptot 


r(10- 12 s) 


r eX pt(10- 12 s) 




2.220 


0.145 


-0.064 


1.027 


3.327 


1.98 


1.23 ±0.09 




2.220 


0.138 


-0.051 


1.027 


3.334 


1.97 






2.220 




-0.110 


1.027 


3.137 


2.10 




n 6 - 


2.217 




-0.127 


1.014 


3.104 


2.12 





suggested in [52,18] that a failure of local duality in the b — > ccs channel, which has smaller energy 
release than that in b — > cud, will further enhance T(B) and suppress Bsl- However, this explanation 
encounters a problem: The charm counting n c will increase and become as large as 1.30 [50], which 
is too large compared to the experimental value n c = 1.12 ± 0.05 [49]. One way out of this difficulty 
for n c is proposed in [53] that a sizeable fraction of b — ► ccs transitions can be seen as charmless 
b — > s processes. In the present paper we will not pursue any of the aforementioned possibilities as 
none of them can explain the lifetime difference between A b and B^- The recipe we are going to 
discuss below [see (3.23)] will solve all the problems with Bsl, F(B), n c and r(Aj,)/r(Bd). 
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It has been advocated in [16] that, unlike the semileptonic inclusive case, since OPE 
cannot be rigorously justified for nonleptonic inclusive decays, the failure of explaining the 
observed lifetime ratio r(Ab)/r(Bd) implies that the assumption of local duality is not correct 
for nonleptonic inclusive widths. It is further suggested in [16] that corrections of order I/tuq 
should be present and this amounts to replacing the heavy quark mass by the mass of the 
decaying hadron in the ??Iq factor in front of all nonleptonic widths. In the following we shall 
see that the ansatz 



nl — > Ini.(^) (3.23) 



will not only solve the short A b lifetime problem but also provide the correct absolute decay 
rates for bottom hadrons. 

Employing the hadron masses 

m Bi = 5279.2 ± 1.8 MeV [42], m B - = 5278.9 ± 1.8 MeV [42], 

m Ab = 5621 ± 5 MeV [54], (3.24) 

we obtain 

r tot (A 6 ) = [3.986 + (0.075 + 0.1055)r] x 10~ 13 GeV + r SL (A 6 ), 
F tot (B d ) = [2.908 + (-0.6445! + 0.739e 2 )] x 10- 13 GeV + r SL (B), 
T tot (B-) = [2.907+ (-0.171Si + 1.974£i)] x 10- 13 GeV + T S l(B), (3.25) 

with Tsl(A6) and Tsl(B) being given by (3.3) and (3.7), respectively. Consequently, 

liM = 0.78 - 0.13ei + 0.15^ - (0.015 + 0.021S)r. (3.26) 
T \ B d) 

Comparing this with Eq. (3.19) we see that the main effect of including linear l/m^ corrections 
is to shift the central value of the lifetime ratio from 0.99 to 0.78. Moreover, the experimental 
value r(A 6 ) / r{B d ) = 0.79 ± 0.06 [12] indicates that the remaining contribution — 0.13s i + • • • 
in (3.26) is at most ±6%. It is also evident from (3.25) that the discrepancy between theory 
and experiment for the absolute decay width of B mesons is greatly improved. 

The most dramatic effect due to the ansatz (3.23) occurs in the lifetime pattern of bottom 
baryons. Employing the bottom-baryon masses (2.36), (3.24) and m Sfc = 5803.7 ± 7.1 MeV, 5 



5 We have used the CDF mass of the A& [see (3.24)] to update the mass prediction given in [39]. 
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some large enhancement to various nonleptonic contributions to the decay widths of bottom 
baryons is shown in Table II. We see that the improved A b lifetime is in agreement with 
experiment and the new hierarchy of bottom-baryon lifetimes emerges as 

r(K) > r(S 6 ") > r(S°) > r(n 6 "), (3.27) 

which is drastically different from the previous one: The longest-lived fib among bottom 
baryons in the conventional OPE now becomes shortest-lived. Needless to say, it is of great 
importance to measure the hierarchy of bottom-baryon lifetimes in order to test the ansatz 
(3.23). The branching ratios of semileptonic inclusive decays are calculated from Table II to 
be: 

B{A b -> Xev) = 8.9%, B(E° b -> Xev) = 7.8%, 

B(E b -> Xev) = 8.4%, B{tt b -> Xev) = 6.9%. (3.28) 

Since serious and precise measurements of the hierarchy of lifetimes of bottom baryons 
may not be available in the very near future, 6 it is thus important to carry out more precise 
measurement of the B s lifetime. An application of the prescription (3.23) will modify the 
prediction [10] 

r(B s ) 



to [16] 



r(B d ) ' 
r(B s 



r(B d 



1 ± 0(1%) 



= 0.938 



(3.29) 



(3.30) 



for the average B s lifetime. The current world average is t{B 8 )/t{B^) = 0.98 ± 0.05 [12]. 

Table II. Various contributions to the decay rates (in units of 10 -13 GeV) of bottom baryons. The 
ansatz (3.23) has been applied to enhance the nonleptonic b quark decay and nonspectator effects. 





pdcc 


pann 


pint 


r SL 


ptot 


r(10- 12 s) 


r cxpt (10- 12 s) 




3.986 


0.260 


-0.116 


1.027 


5.157 


1.28 


1.23 ±0.09 


^0 


4.678 


0.290 


-0.107 


1.027 


5.888 


1.12 






4.678 




-0.231 


1.027 


5.474 


1.20 






5.840 




-0.335 


1.014 


6.519 


1.01 





6 The current LEP results for the lifetime of E b are (1.35lg;^18;^)ps by ALEPH [55] and (1.5±g;J± 
0.3)ps by DELPHI [56]. The average is r(S 6 ) = (1.39loi|)ps. Evidently, the uncertainty is still too 
large to have a meaningful test on the prediction (3.27). 
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IV. LIFETIMES OF CHARMED BARYONS 



In Sec. Ill we see that a replacement of the heavy quark mass with the decaying hadron 
mass in the rriq factor in front of nonleptonic widths provides a much better description of 
the lifetimes of the A & baryon and B mesons. It is claimed in [16] that a much better fit to 
the charmed hadron lifetimes is also achieved if T^l for charm decay approximately scales 
with the fifth power of charmed hadron masses, apart from corrections of order 1/m 2 .. We will 
carefully examine the applicability of this recipe in this section. For a theoretical overview of 
charmed baryon lifetimes, the reader is referred to the review of Blok and Shifman [2] . 

We begin with the semileptonic inclusive decay of the D meson: 

rp-^^WM,)^..^-^.,.)^). (4,) 

We find that the experimental values for D + and D° semileptonic widths [42] can be fitted 
by the quark pole mass m c = 1.6 GeV. Taking m s = 170 MeV, we then have x = {m s /m c ) 2 = 
0.0113 and 

I (x, 0, 0) = 0.9166, h(x, 0, 0) = 0.9556, I 2 {x, 0, 0) = 0.9665 (4.2) 

for charm decay. Repeating the same exercise for charmed baryons, we obtain the charmed- 
baryon semileptonic decay rates 

T(A C -> Xev) = T(E C -> Xev) = 1.533 x 10~ 13 GeV, 

T(n c -> Xev) = 1.308 x 10~ 13 GeV, (4.3) 

which are larger than that of the D meson: 

T(D -> Xev) = 1.090 x 10~ 13 GeV. (4.4) 

The prediction (4.3) for the A c baryon is in good agreement with experiment 

T(A C -> Xev) cxpt = (1.438 ± 0.543) x 10~ 13 GeV. (4.5) 

For charmed baryons S c and Q c , there is an additional contribution to the semileptonic width 
coming from the Pauli interference of the s quark [29]. From (2.5) we obtain 

r int (S c - Xev) = Ir^pecr^l-^Xl+x), 

r int (fi c - Xev) = irU lspcc r nc (l - x 2 )(5 + x), (4.6) 
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where we have applied (2.28) for charmed baryon matrix elements, r' = r /|Kd| 2 and 

r = mo 7 rc»M > ^nspec = lOTT — . (4.7) 

We shall see later that, depending on the parameter r, the nonspectator effect in semileptonic 
decay of S c and Q c can be very significant, in particular for the latter. 

We now turn to the nonleptonic inclusive decays of charmed hadrons. It is well known that 
the longer lifetime of D + relative to D° comes mainly from the destructive Pauli interference 
in D + decay [57,3]. However, it is also known that, depending on the parameters B x and 
especially e±, the Pauli interference T™ t (D + ) in analog to r™ t (S~) given by (3.12) can be 
easily overestimated and may even overcome the c quark decay rate r dcc so that the resulting 
nonleptonic width becomes negative ! This certainly does not make sense. It has been 
discussed in great length by Chernyak [25] as how to circumvent the difficulty with the lifetime 
of D + . We shall not address this issue in the present work and instead focus on the lifetimes 
of charmed baryons. Our purpose is to apply the ansatz similar to (3.23) and see if a better 
description of charmed baryon lifetimes can be achieved. 

In addition to the destructive Pauli interference T 1 "*, there exists another Pauli interference 
term PJ?* in charmed baryon decay which arises from the constructive interference between 
the s quark produced in the c quark decay and the spectator s quark in the charmed baryon. 
Since the expressions of r ann and T 1 ! 1 * for charmed baryons are similar to (3.15) for bottom 
baryon decays, here we will only write down the expressions for T l + described by the last term 
in Eq. (2.4): 



Tf(~ c ) = -Ir or/nspec rH c (l - x 2 ){\ + x)(Bc\ - 2 Cl c 2 - N c c{), 
Ff(n c ) = -ir r/ nS pecrn c (l - x 2 )(5 + x)[Bc\ - 2c lC2 - N c cj). (4.8) 



4 

i] 
6 

It is easily seen that (4.8) is reduced to (4.6) when c± — 1, c 2 = 0, N c = 1 and V ud = 1. The 



and Q c baryons also receive contributions from Cabibbo-suppressed W^-exchange: 
r ann (S+) = |K s /Kd| 2 r ?7nspccr Sc (l - x 2 )(b{c{ + c\) - 2 Cl c 2 ), 



r ann (fi c ) = Q\V us /V ud \ 2 T oVnspcc r nc (l - x 2 )(B(c\ + cl) - 2c lC2 ). (4.9) 
The Q c matrix element [see Eq. (2.18)] 

(n c \(cs)(sc)\n c ) = -6|^;(0)| 2 (2m Cc ) (4.10) 
25 



accounts for the factor 6 in Eq. (4.9). 

To proceed we employ the Wilson coefficients 

din) = 1.35, c 2 (//) = -0.64 (4.11) 

evaluated at the scale /i = 1.25 GeV. From Eqs. (3.21) and (3.22) we obtain B(fi) ~ 
0.74i?(/! had ) ~ 0.74 and r(/i) ~ 1.36 r(/z had ). Repeating the same exercise as the bottom 
baryon case, the results of calculations are exhibited in Table III. We see that the lifetime 
pattern 

r£f) > r(A+) > r(~°) > r(Q°) (4.12) 

is in accordance with experiment. It is evident that when nonspectator effects in semileptonic 
decay are included, as shown in parentheses in Table III, the discrepancy between theory and 
experiment is improved. This lifetime hierarchy (4.12) is qualitatively understandable. The 
S+ baryon is longest-lived among charmed baryons because of the smallness of W^-exchange 
and partial cancellation between constructive and destructive Pauli interferences, while Q c is 
shortest-lived due to the presence of two s quarks in the fl c that renders the contribution of 
T™* largely enhanced. It is also clear from Table III that, although the qualitative feature of 
the lifetime pattern is comprehensive, the quantitative estimates of charmed baryon lifetimes 
and their ratios are still rather poor. 



Table III. Various contributions to the decay rates (in units of 10 -12 GeV) of charmed baryons. 
When nonspectator effects in semileptonic decay are included, the predictions are shown in paren- 
theses. Experimental values are taken from [42]. 





pdec 


pann 


pint 


ri? r SL 


■ptot 


r(10- 13 s) 


r expt (10- 13 s) 


5 


0.903 


0.858 


-0.238 


0.306 


1.829 


3.60 


2.06 ±0.12 


' — 'c 


0.903 


0.042 


-0.226 


0.423 0.306(0.498) 


1.447(1.639) 


4.55(4.02) 


o cr+0.7 
O.O_ 04 


^0 

' — 'c 


0.903 


0.817 




0.423 0.306(0.498) 


2.448(2.640) 


2.69(2.49) 


0.98i°;| 




0.968 


0.224 




1.256 0.262(0.772) 


2.710(3.220) 


2.43(2.04) 


0.64 ±0.20 



In order to have a better quantitative description of nonleptonic inclusive decays of charmed 
baryons, we shall follow [16] to assume that T NL scales with m 5 Hc instead of m^: 



r NL (A c ) : r NL (~+) : r NL (s°) : r NL (n c ) 
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where T^l is the nonleptonic decay rate calculated in the framework of the heavy quark 



expansion and it has the form 

a + b/m 2 c + c/m 3 c + 0(l/m 4 c j\. (4.14) 



(o) _ G 2 F ml 
NL 192tt 3 



To compute the absolute decay width, we introduce a parameter A so that 

r ™_r NL =Ar<»>(^) 5 . (4.15) 

Unlike the ansatz (3.23) for bottom hadrons, it will become clear shortly that A is much less 
than unity for charmed hadrons. Applying the prescription (4.15), treating A, r, B as free 
parameters and fitting them to the data of charmed baryon lifetimes [42], we find 

A = 0.18, r = 1.72, S = 1.46, (4.16) 

where r and B are renormalized at fi = 1.25 GeV. The numerical results are summarized in 
Table IV. Contrary to the previous case, a prefect agreement with experiment will be achieved 
if nonspectator effects in semileptonic decay are not included. 

Table IV. Same as Table III except that the ansatz (4.15) has been applied to enhance the 
nonleptonic c quark decay and nonspectator effects. 



r dec pann pint pint p^ ptot r(lCT 13 s) T cxpt ( l(T 13 s) 



' — 'c 



0.960 2.753 -0.884 0.306 3.136 2.10 2.06 ±0.12 

1.404 0.195 -1.231 1.227 0.306(0.837) 1.902(2.432) 3.46(2.70) 3.5±g;I 

1.415 3.868 1.238 0.306(0.837) 6.828(7.358) 0.96(0.89) 0.98±g;?i 

2.389 1.668 5.775 0.262(1.675) 10.09(11.51) 0.65(0.57) 0.64 ±0.20 



Let us examine the fitted parameters (4.16) in more detail. The value r = 1.72 is fairly 
reasonable as it implies |^ c (0)| 2 = 1.1 x 10 _2 GeV 3 , which is consistent with those of hyperons 
and bottom baryons. Then, does it mean that our previous estimate of r for charmed baryons 
[see (2.34)] is too small ? In our opinion, the enhancement of |^ c (0)| 2 is likely due to the fact 
that \ip%(0) | 2 is not simply equal to ^/|>^d with y — 1 and f D 200 MeV. We conjecture 
that a more realistic value of y is probably close to 3 for charmed baryons and to unity for 
bottom baryons. 

As for the parameter B(/i), it is expected to be less than unity if the valence-quark ap- 
proximation is believed to be valid at a lower hadronic scale. Therefore, it is not clear to us 
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why B(/i) is larger than unity and what is its implication. The smallness of A is attributed to 
the fact that the inclusive nonleptonic decays of charmed baryons are not dominated by the 
c quark decay. Nonspectator effects of W-exchange and Pauli interference terms are expected 
to be of order 



where the factor 16n 2 is a two-body phase-space enhancement relative to the three-body 
phase space of heavy quark decay. Realistic calculations (see Tables III and IV) indicate 
that nonspectator contributions are comparable to and even dominate over the c quark decay 
mechanism. This implies that the charmed quark is not heavy enough (i.e. the energy release 
is not sufficiently large) to make a sensible and meaningful heavy quark expansion. For bottom 
hadrons, we see in Sec. Ill that at least for the A& baryon and B mesons, the nonleptonic 
decay rate is approximated by 



where V is the heavy quark decay rate. However, we find for charmed baryons that 
V dcc (m H Jm c f are 5.36, 7.84, 7.84, 13.24 (in units of 10~ 12 GeV), respectively, for A c , 5+ S° 
and Vt C) where r dec is taken from Table III. Therefore, even in the absence of 1/m 2 and \jm\ 
corrections or even the heavy quark expansion converges, the scaled nonleptonic c-quark decay 
rate T dec (mH c /m c ) 5 already exceeds the experimental decay widths: 3.20, 1.88, 7.72, 10.28 (in 
units of 10" 12 GeV) [42]; that is, 



except for the S°. The presence of large nonspectator contributions (see Tables III and IV) 
will make the discrepancy between theory and experiment for decay widths even much worse. 
Hence, we have to introduce a parameter A 1 to suppress the absolute rates. However, since 
A is an entirely unknown parameter in theory, the recipe of scaling Tnl with the fifth power of 
charmed hadron mass is ad hoc and does not have the predictive power for the absolute decay 
widths. We conclude that, although the ansatz (4.13) provides a much better description of 
lifetime ratios for charmed baryons (apart from the annoying parameter B), the prescription 
(4.15) appears unnatural and unpredictive for describing the absolute inclusive decay rates of 
charmed baryons due to the presence of the unknown parameter A. Since the heavy quark 



167r 2 (A Q c D /m c ) 3 -0.5-0.7, 



(4.17) 




(4.18) 




(4.19) 
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expansion converges very badly, local duality is thus not testable in inclusive nonleptonic 
charm decay. 

V. DISCUSSIONS AND CONCLUSIONS 

We have analyzed the lifetimes of bottom and charmed hadrons within the framework 
of the heavy quark expansion. Especial attention is paid to the nonperturbative parameter 
^baryon four-quark matrix elements for baryons. We found that the large-iV c relation 
A™ eson ~ iV c A2 aryon is satisfactorily obeyed by bottom hadrons. We have followed [10] to 
parametrize the four-quark matrix elements in a model-independent way. Baryon matrix 
elements are evaluated using the NQM and the bag model. The bag-model estimate for 
bottom-baryon matrix elements is smaller than that of the NQM by a factor of ~3. The 
hadronic parameter r defined in Eq. (2.27) is estimated in the NQM to be in the range 0.53 
to 0.61 for both bottom and charmed baryons. Nonspectator effects in inclusive nonleptonic 
decays are then studied in detail. The main results of our analysis are as follows. 

1. Using the charmed quark pole mass fixed from the measured semileptonic decay 
widths of D + and D°, we have calculated I/tuq nonperturbative corrections to 
the semileptonic inclusive widths for other heavy hadrons. We found that while 
TsiX-B) is very close to r S L(Ab), T SL (D) is smaller than r S L(A c ). The predicted 
semileptonic decay rates for the B meson and the A c baryon are in good agreement 
with experiment. This implies that global duality is valid for inclusive semileptonic 
decay. For charmed baryons S c and Q c , there is an additional contribution to the 
semileptonic width coming from the constructive Pauli interference of the s quark. 
This interference effect is sizeable for the S c and becomes overwhelming for the Q c . 

2. The lifetime pattern of the bottom baryons is predicted to be r(f2ft) ~ t(E^) > 
r(Afe) ~ t(2°). Nonspectator effects due to lU-exchange and destructive Pauli in- 
terference account for their lifetime differences. The model-independent expression 
in the OPE for r(A b )/r(_B (i ) is given by (3.19), which is difficult to accommodate 
the data without invoking unnaturally too large values of hadronic parameters. 
Irrespective of the short A b lifetime problem, the calculated absolute decay width 
of the charged B~ meson is at least 20% too small compared to experiment. Since 
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the predicted Tsl(B) agrees with data, the deficit of the B meson decay rate is 
blamed on the nonleptonic width. 

3. Unlike the semileptonic decays, the heavy quark expansion in inclusive nonlep- 
tonic decay cannot be justified by analytic continuation into the complex plane 
and local duality has to be assumed in order to apply the OPE directly in the 
physical region. The shorter lifetime of the A b relative to that of the B d meson sug- 
gests a significant violation of quark-hadron local duality. The simple ansatz that 
Tnl - ► TniX^^/?^) 5 not only solves the lifetime ratio problem but also provides 
the correct absolute decay widths for the baryon and the B meson. The hier- 
archy of bottom baryon lifetimes is modified to 7"(A&) > t(S^) > r(S°) > r(Q&): 
The longest-lived f2& among bottom baryons in the OPE approach now becomes 
shortest-lived. This ansatz can be tested by measuring the E& lifetime in the near 
future. More precise measurement of the B s lifetime provides another quick and 
direct test of local duality. 

4. The lifetime hierarchy r(S+) > r(A c ) > r(S°) > t(Q c ) is qualitatively understand- 
able in the OPE approach but not quantitatively. Apart from an annoying feature 
with the parameter B, a better description of inclusive decays of charmed baryons 
is achieved by scaling Tnl with m b Hc instead of m? c . Contrary to the bottom case, 
a small parameter A 1 has to be introduced, namely Tnl — ► ArNiX^/fc/^c) 5 , 
otherwise absolute decay widths of charmed baryons will be largely overestimated. 
Since A is an entirely unknown parameter in theory, it renders the above prescrip- 
tion unnatural and less predictive. As the heavy quark expansion in charm decay 
converges very badly, it is meaningless to test local duality in nonleptonic inclusive 
decay of charmed hadrons. 

We conclude that the recipe of allowing the presence of linear 1/mg corrections by scaling 
the nonleptonic decay widths with the fifth power of the hadron mass is operative in the 
bottom family but becomes unnatural in charm decay. Can this prescription be justified in a 
more fundamental way ? It is interesting to note that a PQCD-based factorization formulism 
has been developed for inclusive semileptonic B meson decay [23]. This approach is formulated 
directly in terms of meson-level kinematics. Quark-hadron duality can be tested by comparing 
results obtained from quark-level kinematics and those from meson kinematics. The validity 
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of global duality has been demonstrated in the general kinematic region up to 0{\jmq)\ 1/rriQ 
corrections to inclusive semileptonic widths are indeed nontrivially canceled out. When this 
factorization approach is generalized to nonleptonic decays and to heavy baryons, it is natural 
to expect that T^(B) /T^(Ab) ~ (mg /mA t ) 5 if local duality is violated. Since the application 
of PQCD and hence the factorization scheme of [23] to charm decay is very marginal due to the 
fact that the charmed hadron scale is not sufficiently large, the scaling behavior of T NL with 
m H Q occurred in the bottom decay is no longer anticipated in inclusive nonleptonic decays of 
charmed hadrons. 
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